Entering the "magnetic highway": energetic particle anisotropies 

at the heUospheric boundary 



V. Florinski 

Department of Physics and Center for Space Plasma and Aeronomic Research, University 

of Alabama, Huntsville, AL 35899 

Received ; accepted 



- 2 - 



ABSTRACT 



In August of 2012 the Voyager 1 space probe entered a distinctly new region 
of space characterized by a virtual absence of heliospheric energetic charged par- 
ticles and magnetic fluctuations, dubbed a "magnetic highway". Prior to their 
disappearance, the particle distributions strongly peaked at a 90° pitch angle im- 
plying a faster particle escape along the magnetic field lines. We investigate the 
process of particle crossing from the heliosheath region into the "magnetic high- 
way" region using a kinetic approach resolving scales of the particle's cyclotron 
radius and smaller. We show that a "loss-cone" type distribution naturally arises 
as the orbiting particles enter a region of space with an extremely low pitch- angle 
scattering rate. 

Subject headings: cosmic rays — magnetic fields — turbulence — solar wind 
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Introduction 



Since the end of 2004, the Voyager 1 space probe has been exploring the region of 
space downstrearn . of the solar wind t ermination shock, known as the heliosheath (e.g.. 



Stone et al. 



2005 



Decker et al 



is generally in a turbulent state (IBurlaga fc Ness 



20051). Observations sh owed that the heliosheath plasma 



2010l ) with the mean magnetic field 



magnitude increasing with heliocentric distance from some 1 just beyond the shock 
(94 AU) to about 3 /i G by mid-2012 when the spacecraft was at 122 AU from the Sun 



( IBurlaga &: Ness 



2OI2I ). The heliosheath is permeated by energetic ions with energies below 



several MeV, produced inside the heliosph ere from interstellar material (neutral ato ms) and 



accelerated at the termination shock (e.g.. 



Decker et al 



2008 



Florinski et al. 



20091) 



Starting in August of 2012, the HEP intensity ex perienced a s eries of rapic 



and increases after which they disappeared altogether (jStone 



2012 



Decker et al. 



drops 



2012bl ). 



The last event was accompanied by a jump in the magneti c field strength from 3 to 



4 yuG without a change in direction (IBurlaga fc Ness 



2OI2I ). Remarkably, the magnetic 



field became very smooth after crossing the boundary with turbulent fluctuations virtually 
absent. Such a change in the properties of the magnetic fleld was expected to happen across 
the hehopause since it was generally thought that interstellar turbulence existed on vastly 



greater scales (lArmstrong et al. 



I995I ). The absence of magnetic fleld rotation casts some 
doubt on the interpretation that the boundary just crossed is indeed the heliopause. It is 
clear, however, that particle escape along the fleld lines in the new region is very efficient, 
which prompted the name "magnetic highway" (MH). It is not clear at this time whether 
the new region is part of the heliosphere or the interstellar cloud around it, because of 
the nonfunctional plasma instrument on Voyager 1. In the scenario discussed below this 
difference is of no significance as long as the magnetic field remains turbulence free beyond 
the boundary. 
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In this Letter we simulate the process of energetic particle crossing from a turbulent 
region (the heliosheath) into a laminar region (the MH or LISM). Specifically, we propose 
an explanation for the observational result that the intensity of particles streaming along 
the magnetic field decreased faster than tho se gyrating at nearly 9 0° angles resulting in 



a double loss cone pitch-angle distribution (ICummings et al. 



20121 ). Because the HEP 



intensity changes occurred on short time scales (i.e., comparable to VVi/^^, where Vyi is 
Voyager 1 speed and Tg is the cyclotron radius of a typical HEP ion), a model must be 
capable of resolving scales of and smaller. A diffusive approach is obviously invalid here 
because the particle populations were highly anisotropic immediately beyond the boundary. 
Using a fully kinetic approach in a simple numerical model we show that a flattened 
pitch-angle distribution is a natural consequence of gyrating particles traveling farther from 
the boundary before escaping under scatter-free conditions. 



2. The kinetic transport model 

The large-scale geometry of our model is illustrated schematically in Figure [H We 
consider a narrow region immediately adjacent to the boundary on either side. The x-axis 
is normal to the boundary and the z-axis is in the direction of the mean magnetic field. It 
is assumed that the magnetic field lines in the MH region are pressed against the boundary 
for a distance Zmax (the "connection region" in Fig. [1]) after which they separate from the 
boundary. A particle reaching the end of the connection region is assumed to escape the 
region and is removed from the system. 



The backgrou nd plasma velocity indirectly measured by Voyager 1 was very small 



during 2010-2012 (IKrimigis et al. 



2011 



Decher et al. 



2012al ). For simplicity we assume 



here that the plasma velocity is zero on both sides of the boundary (this assumption is 
not critical as long as the plasma speed is much less than the particle velocity f , which is 




Fig. 1. — An illustration of the particle transport model geometry. We consider a narrow 
"connection region" where the magnetic field presses against the boundary. 
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usually satisfied everywhere in the heliosheath) . We take the mean field to be -Bhs = 3 fiG 
and -Bmh = 4 fiG on the heliosheath side, and the MH side of the boundary, respectively. 
The general Boltzmann equation describing energetic particle transport in a magnetized 
plasma with turbulent fiuctuations may be written as 



df 



df , df 



- + . v/1 - /^^ cos ^- + v^^- -n- = Sf, 



where /(x, z, /i, y?) is the phase space density, /i is the pitch-angle cosine, ip is the gyrophase, 
and fl is the cyclotron frequency. The operator S, describing particle scat tering in solid 



Jokipii 



1972|) 



angle, could be obtained from the theory of wave-particle interaction (e.g. 
The transport across the magnetic field is produced entirely by scattering of the particles 
in gyrophase. Because the plasma is stationary, there is no electric field and the energy of 
the particles is assumed to be conserved. 



For simplicity we adopt a model with an isotropic scattering operator. 



d 



'1 2N^/ 



+ 



1 — /i^ d(p'^ ' 



(2) 



where D is the scattering coefficient. The parallel mean free path corresponding to 
(2) is A|| = fr, where r = (2D)~^ is the scattering rate. For large r2r used here the 



transport is predominantly along the magnetic field (e.g.. 



Forman et al 



19741 ). Equation 



1) is solved numerically a . 



(e.g. 



MacKinnon fc Craig 



■ong 



t he characteristics using a st ochastic integration method 



1991 



Florinski fc Pogorelov 



2009). Instead of 



u sing E g. (2) 



Ellison et al. 



fll990l ) that 



directly, we use the random walk on a sphere technique of 
avoids the singularity at yU = ±1. For small scattering angles the process is identical to a 
two-dimensional random walk on a plane. 
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3. Intensity and pitch-angle distributions 

We performed a series of simulations varying the scattering rates, D^s and Z^mh, and 
the length of the connection regions Zmax? within a sufficiently wide range. Figure [2] shows 
a typical result for 1 MeV protons obtained with Dhs = 10^'^ s^^, -Dmh = 2 x 10^^ s^^, and 
-^max = 32 AU. Such a large difference in scattering rates is required to produce a noticeable 
difference between particles traveling at different pitch angles. The MH region is essentially 
scatter-free (Ay = 23 AU) and particles are only able to penetrate a distance comparable to 
Vg = 0.0024 AU before escaping into the LISM. Intensity profiles of yU = 1 and fi = ions 
are shown with a solid and dashed lines, respectively. 

The particle distribution is assumed to be isotropic at the left (heliosheath) boundary 
of the simulation and the two curves are almost identical up to x = 0. Afterwards, ions 
streaming along the magnetic field diminish earlier than ions near 90° pitch angle. The right 
(MH) boundary plays no role in the simulation because all particles exit the simulation box 
before reaching it. The difference in penetration depth is of the order of a cyclotron radius 
in the MH. We note that the result remains essentially unchanged if the ratio -DMn/^max is 
kept a constant. This means that the model applies for larger or smaller scattering rates as 
long as the unknown escape distance is adjusted accordingly. 

Figure [3] plots the HEP pitch-angle distributions on the heliosheath side at x = —0.005 
AU (solid line) and on the MH side at x = 0.005 AU (dashed line). The initially isotropic 
ion population becomes highly anisotropic past the boundary as it becomes depleted in 
/i = ±1 particles. Some fi = particles have also scattered and escaped from the system, 
so that the total intensity is significantly reduced. This result is consistent with Voyager 
1 observati ons showing a fiattened (pancake like) pitc h-angle ion distribution beyond the 



boundary (ICummings et al. 



2012 



Decker et al. 



2012bh . 
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Distance from the boundary, AU 



Fig. 2. — Intensity profiles of /i — 1 (solid line) and /i — (dashed line) particles as a 
function of distance from the boundary (the dotted line). 
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Fig. 3. — Pitch-angle distributions of HEPs in the hehosheath {x — —0.005 AU, sohd hne) 
and in the MH region (x — 0.005 AU, dashed hne) The intensity is normahzed to one at the 
left boundary of the simulation. 



- 10 - 



4. Discussion 

In this Letter we proposed a simple model to explain the double loss cone type 
distributions of MeV heliosheath particles observed immediately beyond the boundary 
crossed by Voyager 1 in 2012. The model is based on a concept that particle scattering 
rates are much smaller in the MH region than in the heliosheath. Particles gyrating near 
fi = 1 are able to penetrate deeper into the new region because of their larger gyro-radius, 
whereas particles streaming along the field rapidly escape into the LISM. The model does 
not require magnetic mirrors on either side of the Voyager crossing point, but only that 
there exist a boundary magnetic connection region of co nsiderable extent (a few t ens o " 



AU). Our model is unrelated to the concept proposed by lMcComas fc Schwadrod (120121 ). 
where the drop in energetic particle flux is a consequence of Voyager 1 leaving the region 
magnetically connected to the termination shock. 

Our very simple model qualitatively reproduces the observed changes in particle 
distribution across the boundary. The same model may be applied to ions with higher 
and lower energies by adjusting the scattering rate accordingly. The contrast between 
particles traveling at and 90° pitch angles in our model appears to be less than observed. 
This is primarily due to a strictly perpendicular magnetic field geometry used here. The 
model therefore lacks diffusive transport across the magnetic field from magnetic field 
line meandering. A proper statistical treatment of this effect on kinetic scales requires a 
prescription of the complete three-dimensional turbulent magnetic field, which is beyond 
the scope of this paper. Magnetic field line diffusion would make particle intensity decrease 
at a slower rate. Consequently, the HEPs may be still around several Vg deep into the MH 
region. 

Another simplification is our use of a very thin boundary between the two regions. 
The boundary could be diffuse and/or turbulent, which would also make the HEP intensity 
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decrease slower than predicted here. It should be pointed out that the true radial extent of 
the particle decrease region is unknown because of a large uncertainty in the background 
plasma speed measurement. If the flow was slightly anti-sunward, the region could be very 
thin (perhaps of the order of r^). Irrespective to this result, we think that the basic premise 
of our model (particle escape into the LISM at different rates) adequately describes the 
fundamental physics of the newly discovered boundary. 

This work was supported, in part, by NASA grants NNX10AE46G, NNX12AH44G, 
NSF grant AGS-0955700 and by a cooperative agreement with NASA Marshall Space Flight 
Center. 
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